Dynamics of the Ising Spin-Glass Model in a Transverse Field 
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We use quantum Monte Carlo methods and various analytic approximations to solve the Ising spin- 
glass model in a transverse field in the disordered phase. We focus on the behavior of the frequency 
dependent susceptibility of the system above and below the critical field. We establish the existence 
of an exact equivalence between this problem and the single-impurity Kondo model at the quantum 
critical point. Our predictions for the long-time dynamics of the model are in good agreement with 
experimental results on LiHoo.i67Yo.g33F4. 
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The physics of frustrated quantum spin systems is a 
fascinating and rapidly growing area of condensed matter 
physics . One of the most widely investigated systems 
is the Ising spin glass in a transverse field, a model that 
combines simplicity and experimental accessibility 
The Hamiltonian of the model reads, 



H 



(1) 



where , ii — x, z are components of a three dimensional 
spin-1/2 operator at the i— th site of a fully-connected 
lattice of size N and the first sum runs over all pairs of 
sites. The exchange interactions Jij are independent ran- 
dom variables with a gaussian distribution of zero mean 

1 /2 

and variance J = (^Jf^^ , and F is an applied magnetic 
field transverse to the easy-axis z. For F = Eq. |l| 
is the classical Sherrington-Kirkpatrick spin-glass model 
that has a second order phase transition at Tg = J/4. 
When F is finite quantum fluctuations compete with the 
tendency of the system to develop spin-glass order. As 
a result a boundary F(r) appears in the F — T plane 
between spin-glass (SG) and paramagnetic (PM) phases. 

The model of Eq.|l] is relevant for the compound 
LiHoo.i67Yo.833F4, a site-diluted derivative of the 
dipolar-coupled Ising ferromagnet LiIIoF4 An exter- 
nal magnetic field perpendicular to the easy axis splits 
the doubly degenerate ground state of the Ho'^"'' ion. 
This splitting is proportional to and plays the role 
of F in Eq. | |,|. Experimentally, LiHoo.i67Yo.833F4 
is paramagnetic at all temperatures above a critical field 
« 12 kOe [|,|. Below and for T > 25 mK the 
behavior of the non-linear susceptibility indicates a sec- 
ond order transition line between SG and PM phases 
that ends at r° = 135 mK and Ht ^ ^. Investiga- 
tion of the long-time dynamics of this system above this 
line has revealed the existence of a fast crossover in the 
field dependence of the absorption at very low frequencies 
|,|. This crossover is characterized by a steep increase 
of x"('^ 0) across an almost flat line in the F — T plane 
at F «i Fc and up to T - T°. 



While the phase diagram of the model has been stud- 
ied theoretically using a variety of methods j|,^,|l^ , much 
less is known about its dynamics that has only been dis- 
cussed for large F/J |]ll]| , and near the quantum critical 
point 

In this paper we use a recently developed quantum 
Monte Carlo method (QMC) jl^] to find the exact nu- 
merical paramagnetic solutions of the model throughout 
the F — r plane, and also obtain analytic expressions 
that we derive in various limiting cases. We find that the 
behavior of x"{^) above and below the critical field is 
qualitatively different. For F > Fc the zero-temperature 
spectrum of magnetic excitations has a gap A |l^ that 
vanishes as F ^ Fc « 0.76J. At finite but low tem- 
peratures, T < A, the gap edge develops a tail of ex- 
ponentially small weight. On the other hand, for small 
F and low T, we find a narrow feature around w = 
whose intensity decreases rapidly with increasing field or 
temperature as spectral weight is transferred to higher 
frequencies. We further demonstrate that at the quan- 
tum critical point the problem can be exactly mapped to 
the single-impurity Kondo model. The low-energy prop- 
erties of the system in the neighborhood of this point are 
characterized by a new energy-scale. To « 0.08J. At fi- 
nite temperature there is a crossover between the regimes 
just mentioned that is essentially controlled by F up to 
T ^ T°. Wc finally present detailed predictions for the 
F and T-dependence of the low-frequency response that 
are in good agreement with the experimental results on 

LiIIOo.l67Yo.833F4 ||] • 

Bray and Moore [|l^ have shown that the quan- 
tum spin-glass problem can be exactly transformed 
into a single-spin problem with a time-dependent self- 
interaction Q{t) determined by the feedback effects of 
its coupling to the rest of the spins. As we have shown 
elsewhere for a related problem much progress can 
be made by eliminating the self-interaction in favor of an 
auxiliary fluctuating time-dependent field 7y(T) coupled 
to the spins. The free-energy per site J- in the paramag- 
netic phase can then be written as 
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(iT = min I ^ f f drdr'Q^T - r') - In Zi^dQ] \ , 

Q{r) [ 4 Jo Jo J 

Zioc[Q]= J 2?r;cxp -^J^^ jyTdT'Q-\T,T')rjiT)rjir') 



xTrT exp 



dT [j77(r)5"(T) +r5"(T)] 



, (2) 



where T is the time-ordering operator along the 
imaginary-time axis < t < /3. Zi^c can be thought of 
as the average partition function of a spin in an effective 
magnetic field hcs = J'ri(T)ez + re^; whose z-component 
is a random gaussian function with variance Qij). The 
latter is determined by functional minimization of (|^) 
which gives the self-consistency condition |14 



Q{r) 



(3) 



where the average is taken with respect to the probabil- 
ity density associated to Z\oc- We solved Eqs. |^]3|itera- 
tively using the QMC technique that we have described 
elsewhere |13 . The imaginary time axis is discretized 
in up to L = 128 time slices with Ar = (3/L < 0.5. 
An iteration consists of at least 20,000 QMC steps per 
time-slice and self-consistency is generally attained af- 
ter about eight iterations except very close to the quan- 
tum critical point. We mapped the spin-glass transition 
line in the T — T plane using the well known stability 
criterion |14| 1 — Jxioc where xioc = dTQ{T) is the 
local spin-susceptibility. We found a second-order tran- 
sition line ending at a quantum critical point at T = 
in agreement with previous work P,p|pO|] . Going down in 
temperature to T ~ 10"^ J and extrapolating the results 
to T = we determined a precise value for the critical 
field Fc/J = 0.76 ± 0.01, which lies in between previous 
estimates [|2|jl^. 

We have studied the dynamical properties of the para- 
magnetic state throughout the T — T plane, even below 
Tg where it is unstable. Indeed, the analysis of the evo- 
lution of the paramagnetic solution for small T provides 
insight on the physics of this problem as the sates below 
and above Tg are continuously connected. In Fig.^ we 
show the correlation function for several values of F and 
T. For F > Fc (panel a), Q{t) decays exponentially with 
a time-constant tq ~ F^^ ^ (3 that depends only weakly 
on T. This behavior is characteristic of the existence of 
a gap A ~ F in the excitation spectrum of the system. 
For F < Fc and T <^ J (panel c), Q{t) also decays very 
rapidly for short times, r ^ tq ~ /3~^. For r ^ tq, 
however, it exhibits a very slow variation which indicates 
that the presence of excitations in the low-energy end of 
the spectrum, w <C T. With increasing temperature tq 
increases and reaches a value 0{J~^) when T ^ J. It 
then becomes no longer possible to distinguish two dif- 
ferent time scales. The case F = Fc at low temperature 
(panel b) is intermediate between the other two and the 



long-time behavior becomes a power-law, Q{t) cx as 
T ^ 0, as anticipated by Miller and Huse using in- 
ternal consistency arguments. The solid lines in Fig. |l| 
are the results of various analytic approximations that 
we discuss next. We begin by considering T / J ^ 1. In 
this case, the effective fields appearing in Eq. ^ are dom- 
inated by their a;— component. We may thus evaluate the 
trace of the time-ordered exponential under the integral 
in Z\oc using a low-order cumulant expansion. To the 
lowest non-trivial order we find. 



(4) 



where rrix = l/2tanh(/3F/2) is the zero-th order trans- 
verse magnetization. Eq. ^ predicts the existence of a 
gap A = VF^ — 2F Jma; in the excitation spectrum. This 
gap has a very weak temperature dependence for T <C F 
and, at T = 0, vanishes at Fc/J = 1, overestimating 
the critical field. It can be shown ||l^ that, for T ^ A, 
the gap-edge develops a tail that carries a small weight 
0{e-^'^). 
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FIG. 1. Q(r) as function of r//3. The crosses correspond 
to QMC data. The error bars are smaller than the size of 
the crosses. The solid lines are obtained using the analytic 
expressions discussed in the text. 

This procedure breaks down sXT / J ^ 1 and a different 
approach must be taken in order to describe the physics 
at low fields. At F = 0, the problem reduces to the 
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classical Sherrington-Kirkpatrick model and Q{t) = 1/4 
at all temperatures. We thus expect that for F/J ^ 1 
the effective field /loff(wn) will be dominated by its cj = 
component. Indeed, setting for the moment r]{Lu) = for 
w ^ in Eq. ^ and performing the functional average 
solely over static fields '^{t) =770, we find 



Q{uj„ 



F2 + .Pril 



Tm{rja) cos 



(5) 



where m{rjQ) — 1/2 tanh(/3-\/T2 + J2?yJ/2) and cosO^ = 
r / ^/r^+J^rj^ . The average is performed with 
respect to the probability distribution, V{rio) oc 
exp(-/37;2/2xioc) cosh(/3Vr2 + J^Vo/^)- Using Eqs. | 
and |5| we can deduce the region of validity of this ansatz 
from the estimates \r]{ujn)/r]{0)\ - 32FTVj3 for T < J, 
and ^ F/(27r?T,T) when T ^ J. Within this approxima- 
tion, the imaginary part of the response on the real axis 
is given by. 



F2 + J^7% 



5{lo) 



no 



F2 tanh(/3|cj|/2)7'(Vt^2 _ r2) 



2cj2 



F2 
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In this regime the relaxation function x" splits into 
two contributions, an elastic peak at w = and a con- 
tinuum starting at w = F. The fraction of spectral 
weight contained in each of these two contributions is 
determined by the transverse field and the temperature. 
Evaluating the coefRcient of the (5-function in Eq. ^ we 
find that the relative intensity of the central peak varies 
between 1 - (8Fr/j2)2 for T < J, and J2/(4F2 -|- J2) 
for T ^ J. The low-energy states represented by the 
(5-function are responsible for the slow decay observed in 
the long-time behavior of Q{t) at low T in Figj^c. The 
fact that the central peak has zero width is a shortcom- 
ing of the approximation leading to Eq. || as it does not 
capture the slow relaxational processes which broaden it 
Nevertheless, the excellent agreement between the 
numerical and analytic results indicates that the width of 
the central peak must be, in any case, much smaller than 
the temperature. The high-energy states described by 
the inelastic part of the response control the exponential 
decay observed for short times. The decay rate predicted 
by Eq. | is t^^ ~ for T < J, and t^^ - J for 

T ^ J in agreement with the numerical results. 

The approximations discussed above can not be used 
near F = Fc. However, one can still gain insight on the 
dynamics in the critical region by exploiting an interest- 
ing analogy between the model (0) at the quantum criti- 
cal point and the single impurity Kondo problem that we 
establish next. We first perform a Trotter decomposition 
of the time ordered exponential in Eq. ^ and introduce 
intermediate states |cr)((T| at each imaginary time-slice t. 
The trace is now evaluated using the expression 



(7) 

valid in the limit when the witdth of the time slice 
At —>■ 0. The partition function can then be expressed 
in terms of a sum over "histories" , each of them defined 
by a particular sequence of the eigenvalues (J^{t) = ±1/2 
of the intermediate states. Going over to the continuum 
limit and performing the gaussian integral over the aux- 
iliary fields ?7(t), we obtain 



Z\oc = / 2?o-^ exp 



J2 



"'0 

ln(F/2) X (number of spin flips) , 



(8) 



where by "number of spin flips" we mean the number of 
times that the function cr^ (r) changes sign in the interval 
< T < /3. This expression has a form analogous to that 
of the partition function of the single-impurity Kondo 
model in the Anderson- Yuval formalism |16|. In Eq. ||, 
F and J^Q{t) play the roles of the spin- flip coupling J± 
and the long-range Ising-like effective interaction in the 
Kondo model |l^ . Still, there is an important difference 
between the two problems. In the latter the Ising-like in- 
teraction is given, and behaves as (2 — e)/T'^ with e oc 
as T ^ [Q. When the dynamics of the impurity- 
spin is controlled by the strong-coupling fixed point, its 
time-dependent correlation function is also oc t~2 long 
times. In contrast, in our problem, Q{t) is a priori un- 
known. However, one realizes that, if for some F the 
asymptotic behavior of Q{t) is ~ t~2 then, by virtue of 
and the identifications just made, the two problems 
become equivalent at that value of the field and T = 0. 
This is indeed the case at Fc where, as shown in Fig. 
|l|b, the self-consistent solution of the problem at low T 
can be accurately fitted by the finite-temperature gener- 
alization of Anderson and Yuval's Ising-like interaction 
Ip^ using To ~ F~^ as a short-time cutoff. This analogy 
between our problem and the Kondo model provides a 
simple way to estimate the energy-scale associated with 
the low-energy excitations at the quantum critical point. 
It is well known [0 that the local susceptibility of the 
Kondo impurity is given by T^XiociT) — (j){T /T^)), where 
is a universal function with <^(0) « 0.0796 [0 and 
To is the Kondo scale. Since in our case Xioc = J^^ at 
the quantum critical point, it follows that Tq « 0.0796 J. 
We expect the quantum critical region to extend up to 
a temperature Tqc of this order. Our detailed numeri- 
cal results for the temperature dependence of the local 
susceptibility are consistent with Tqc ~ To/4 | p^ . 

Finally, we would like to discuss the experimentally 
observed crossover dynamics in LiHoo.i67Yo.833F4 in the 
light of the results presented here. By measuring the 
response of the system at 1.5 Hz, the experiment [0 is 
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probing the intensity of the low-lying excitations repre- 
sented by the (5-function peak in Eq. ^ 
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FIG. 2. (a) Spectral weight of the central peak in Eq. ^ as a 
function of (F/Fc)^^'^ for several temperatures, (b) Imaginary 
part of the linear susceptibility of LiHoo.i67Yo.g33F4 at 1.5 Hz 
as a function of Ht- Data correspond to T=98, 82, 66 and 25 
mK, from left to right. From Ref. 

In Fig. H we compare the intensity of the latter (upper 
panel) and the experimental data of Ref. [|j (lower panel) 
at various temperatures. To make the connection with 
the experiment we plot the theoretical results as a func- 
tion of (r/Fc)^/^ since the splitting of the ground-state 
doublet of the Ho^+ ion is proportional to the square of 
iJf The temperatures chosen for our calculations are 
such that the ratios T/T^ correspond to the experimen- 
tal values. There are no free parameters other than the 
overall scale of the y-axis. The theoretical curves end at 
the field rc(r) where the condition Jxioc = 1 is fulfilled. 
For lower fields, the system enters the spin-glass phase 
and the experimental data become field-independent. As 
the figure demonstrates the overall behavior of the the- 
oretical curves is in remarkable agreement with the ex- 
periment. Two important features are worth noticing. 
Firstly, the value of the field at the onset of absorption, 
Ht ~ is almost temperature-independent which ex- 
plains the flatness of the experimental crossover line. Sec- 
ondly, for T ^ T° absorption starts well above Tc{T), 



meaning that precursor low-lying excitations appear in 
the paramagnetic phase long before the system freezes. 

In conclusion, we have presented an exact numerical 
solution of the infinite-range Ising spin-glass model in a 
transverse field in the paramagnetic phase. We worked 
out analytic approximations in different limiting cases 
that allow for a physical interpretation of the numeri- 
cal data. In particular, we established for the first time 
an interesting connection between this problem and the 
single-impurity Kondo model. Our prediction for dynam- 
ics of the model are in good agreement with experimental 
results on LiHoo.i67Yo.833F4. 
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